Faddeev-Hopf knots: Dynamics of linked un-knots 
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Abstract 

We have studied numerically Faddeev-Hopf knots, which are defined as those unit-vector 
fields in R 3 that have a nontrivial Hopf charge and minimize Faddeev's Lagrangian. A given 
initial configuration was allowed to relax into a (local) minimum using the first order dissipative 
dynamics corresponding to the steepest descent method. A linked combination of two un-knots 
was seen to relax into different minimum energy configurations depending on their charges and 
their relative handedness and direction. In order to visualize the results we plot certain gauge- 
invariant iso-sur faces. 

The theory of knots has a long history in mathematics and physics and research on this 
topic has again become very active, partially due to the possible real world applications 
(DNA knotting etc). In the late 19th century knots were proposed as models of atoms, 
and although this idea did motivate the mathematical study of knots, it did not last as a 
model for physical reality. Now, with the advent of string theory, the possibility of knotted 
elementary structures has been considered again even in particle physics. In order to give 
knots a physical meaning one has to give them some physical properties, e.g., thickness 
and energy. For example, one may assume that there is a charge distribution on the wire 
and then find its minimum energy configuration 

The starting point for the knots discussed in this paper is classical field theory with 
a unit-vector field denned everywhere in the 3-dimensional space. The unit vector can 
be considered as a point in S 2 ; furthermore, if at spatial infinity all vectors point to the 
same direction (lets say to the north pole, i.e., n = (0, 0, 1)) we can compactify R 3 to S 3 
and the vector field provides a map n : S 3 —>■ S 2 . Such maps can be labeled by their Hopf 



charge, ^(S" 2 ) = Z. As an example we mention the following map with Hopf charge one: 



4(2xz - y(x 2 + y 2 + z 2 - 1)) 
(1 + x 2 + y 2 + z 2 ) 2 ' 
A{2yz + x[x 2 + y 2 + z 2 - 1)) 

(1 + x 2 + y 2 + z 2 ) 2 ' ^ j 

S(x 2 + y 2 ) 
(l + x 2 + y 2 + z 2 ) 2 ' 

(It is easy to see that n — > (0,0, 1) as x 2 + y 2 + z 2 — ► oo, as required). A knot can 
be associated to the vector field by considering how it turns around. In particular, one 
can look at the pre-image of the south pole n = (0,0, —1); normally it is a closed curve 
which may form a knot. For the vector field (|1|) the pre-image of the south pole is the 
ring z = 0, x 2 + y 2 = 1, justifying the name "un-knot". Later it will become clear that 
the curve forming the (un-)knot has also direction and handedness, and the relevant knot 
theory is therefore that of "framed links" 0. [An un-knot of charge Q is obtained from 
(H) if we replace x — > dt[(x + iy)®]/r®~ 1 ,y — > S[(x + iy)Q\/rQ~~ l , where r = \fx 2 + y 2 .} 

In addition to topology we need physics, in particular we need to define an energy 
for the field configuration. One requirement is that it should not be feasible to minimize 
energy by scaling the knotted configuration to zero or infinite size. One early proposal 
was 0] E — ~ J[((9 M n) 2 ] zd 3 x (//-sum over 1 to 3), but this does not contain a scale. In ||] 
Faddeev proposed the Lagrangian 

L = \J {(^n) 2 + gF 2 u } d 3 x, = e abc n a d^n b d u n c , (2) 

as the energy functional. (The topological term can also be written as F 2 V = |(<9 M n x 
d u n) 2 = |[(<9 M n) 2 (<9j,n) 2 — (c^n • d u n) 2 }.) What is important here is that there is now a 
scale defined by the dimensional coupling constant g. Indeed, the derivatives mean that 
the integrated kinetic term scales under r — > Xr as A, and the F 2 term as A^ 1 . Thus 
a given initial knotted field configuration (e.g. ([]])) should evolve (with suitable added 
dynamics) to a stable minimized configuration of a given size and shape. Vakulenko and 
Kapitanskii || obtained the lower limit (VK bound) 

E>c^g-\Q\-*, (3) 

for the energy given by (Q), here c is some constant, g the coupling constant in ([|), and 
Q the Hopf charge. 

The first numerical results with Lagrangian @ were obtained by Faddeev and Niemi 
0. (for popularizations see ||) and Gladikowski and Hellmund 0. They studied mainly 
charge 1 and 2 un-knots, although the possibility of real knots was briefly discussed in 
H 0] . The purpose of this letter is to describe the dynamics of linked un-knots. We find 
that under minimization the linked un-knots can deform to other configurations following 
the rules applicable to framed links (Kirby moves). 



n 1 = 
n 2 = 
n 3 = 
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The results of || [j]] were obtained by a computational method that is different from 
the one used here. Furthermore, in those papers a rotationally symmetric ansatz was used, 
because it reduces the computations to the two-dimensional space containing the z axis. 
However, this reduction also means that stability under rotation-dependent perturbations 
cannot be probed. It turns out that such effects can be seen already for the charge 3 un- 
knot for which the minimum energy configuration is not rotationally symmetric |10| . 
Our method is well suited for knots of more complicated topology and we will report 
some fairly complicated minimization dynamics for linked un-knots. In || |7j another 
simplification was made in handling the unit-vector field. After all, the three components 
of a unit vector are not independent and one can instead use the two components defined 
by w = U + iV = (n 1 + in 2 ) /(l + n 3 ). Unfortunately there is a numerical problem in using 
w instead of n, due to the fact that the range of U and V is from — oo to +00. The singular 
values are attained along a curve forming the knot, i.e., where n = (0,0, —1). It is clear 
that in a relaxation process it is very expensive to try to move infinities. Indeed, in || [7| 
the infinite inertia of the w = 00 ring was compensated by changing the coupling constants 
to best fit the virial theorem. The behavior near this infinity ring is also doubtful, because 
controlling overflow situations may cause unexpected effects. 

To avoid the problems mentioned above we decided from the very beginning to put the 
system on a cubic lattice and use unit-vector fields themselves. Thus, the unknowns were 
n^- fc , where i,j, k give the location on the lattice and a the component, ELiKfc) 2 = 
1, V2, j, k. For the kinetic energy term derivatives of n were discretized on links, that is 



J (d„n(x, y, z)) 2 d 3 r -> £ [(K +ljk - n a ijk ) 2 + (nf j+lk - nf jk ) 2 + (nf jk+1 



n ijkJ 



i,j,k,a 



5. 



For the potential term we discretized as defined in (||) on the v) plaquette. To 
make it completely symmetric we used in F xy , say, 

n ° — y ~^(- n< ijk + n T+ljk + n ?+lj+lk + n ?j+lk)i 

d x n a — > -^[( n f + ijk ~ n ?jk) + ( n ?+ij+ik ~ n ij+ifc)L 
d y n a — > — — n°- fc ) + (nf +1 j +lk — n^ +l - k )}. 

Minimization of the energy was done using the steepest descent method: n new = 
n-oid + dtVL, which corresponds to the dissipative dynamics n = VL. This was sometimes 
accelerated by taking into account also the gradient of the previous step. The expression 
for the gradient was calculated from the discretized Lagrangian using the symbolic algebra 
program Reduce. We did not use Lagrange multipliers, but rather renormalized n new to 
preserve unit length. On the lattice boundary the vector was fixed to n = (0, 0, 1). The 
fact that the system was put into a box with fixed boundaries introduced some boundary 
pressure. For this reason we cannot exactly attain the prediction of the virial theorem 
stating that the topological and kinetic energies should be equal. 
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We did first some calculations on a 50 3 and 100 3 lattices, but it soon became clear 
that even for un-knots the discretization effects could sometimes change the outcome: 
If the lattice is too small some neighboring vectors could turn anti-parallel and break 
the topology. The present computations were done on a 240 3 lattice divided into 4 3 
processors on Cray T3E parallel machine. Each round of iteration took about 2 seconds. 
Typical runs took around 60-70.000 iterations, i.e., 30-40 hours (for charge 1 and 2 linked 
un-knots). Also, to make sure there were no topology-breakup effects we followed the 
nearest-neighbor differences of the vectors, in our computations max|n-n nn — 1| was 
large only at the very beginning and even then the differences never exceeded 0.1 and 
later were typically 0.02. 

It is, in general, hard to visualize vector-fields in three dimensions. Often we associate 
flow lines to vector fields by imagining how a test particle would move if the vector field 
would be its velocity field. In the present case, however, this method is not relevant. The 
point is that there is a global gauge invariance and field lines change in an essential way 
under such global rotations. 

As was said before, the unit vector field provides a map R 3 ~ S 3 —>■ S 2 . We have one 
fixed direction, that of the vector at infinity (n = (0,0, 1)), which defines the north pole 
of S 2 . The unit vector at any other point is defined by its latitude and longitude. In our 
illustrations we have chosen to display the pre-image of some latitude circle, that is, the 
surface where the vector field has a fixed angle with the vacuum direction, i.e., where n 3 
has a fixed value. 

The longitude is described by colors on the iso-latitude surface. This is suitable because 
colors are commonly put on a circle anyway: red — > yellow — > green — > blue — > magenta 
— > redQ The actual position of the meridian is not important, because it is a gauge 



dependent quantity. The figures were made using the program "funcs" |l 1] running on a 
Silicon Graphics 02. 

The vector field (|1|) is not that suitable as an initial configuration for numerical analysis 
because its kinetic energy density decreases only as r -4 as r —>■ oo and therefore the 
integral itself as r _1 . Thus, we introduced some radial squeezing for large r intended 
to minimize interaction with the box boundaries, and some expansion for small r to 
decrease the rapid changes at the center. From such a starting configuration a single un- 
knot relaxed fairly rapidly (in about 30-40.000 iterations) into the minimum configuration. 
Usually we followed the minimization until the changes in the energy were of the order 
10~ 7 per iteration. The same was done for the charge 2 un-knot. In both cases the final 
configuration seemed to be rotationally symmetric, as can be seen from Figure [1|. From 
these results we infer that for our normalizations the constant c in the VK bound (Q) is 
about 380. 

In order to test stability under non-rotational perturbations we took the final result of 



1 Unfortunately the coloring algorithm we used did not completely support the color-circle idea: For 
example, even if the values and 1 were to stand for the same color, a triangle with corner values 0, 
and 1 was painted with the color corresponding to the average 1/3. 
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an un-knot (with charge 1 or 2), cut it into two semi-tori, separated them by the distance 
of 20% of the lattice length by repeating the middle configuration, and embedded the new 
configuration back into the original lattice. The relaxation quickly took this configuration 
into the original circular form Jl2| . 



From Figure |l| one can clearly see that the un-knot has both handedness and direction. 
The handedness can be determined e.g. by holding the ring and checking with which hand 
the thumb is aligned along the color boundaries. A right-handed ring (corresponding to 
(0)) has positive Hopf charge. We can arbitrarily associate a direction to the positive 
charge un-knot, we use the color direction red — > yellow — > green — > blue — > magenta — > 
red. It turns out that for the left-handed (=negative charge) un-knot we should define 
the direction by the opposite color sequence. 

Our main object was to study the relaxation dynamics of linked un-knots. The first 
problem was to construct suitable initial configurations. For a single un-knot we could 
use analytical formulae, but for linked configurations this cannot be as easily done with 
one formula. Besides, it would be better to use the already relaxed numerical results as 
starting configurations. We therefore chose a "cut and paste" technique. 

Two un-knots each in a box of size L 3 (scaled, if necessary, to the same g value) are 
immersed into a box of size (2L) 3 as described in Figure 0. Assume that the ring of the first 
un-knot is on the (x, y)-plane and at infinity the vectors point to the positive z-direction. 
From this single un-knot configuration n we use the part < x < L, < y < L, \L < 
z < \L. The new configuration N in the large cube contains, at first, eight pieces defined 
by n as follows: (the variable range is always < x < ~L, < y < |L, |L < z < |L) 

N(x,y + \L,z + \L) = n(x,y,z), 
N(z, y + \L, z + \L) = n(x,±L,z), 



N(x,y+lL,z + lL) = n(x,y + \L 



z 



N(x+±L,y + lL,z + ±L) = n(|L, y + |L, z), 

N(x + L,y + lL,z + \L) = n(x + \L, y + \L, z), 

N(x + L,?/ + fL,2 + |L) = n(x + lL, \L,z), 

N(x + L,y + \L, z + \L) = n(x + \L, y, z), 

N(x+\L )y + \L,z + \L) = n(±L,y,z). 

The other un-knot is inserted similarly, but first it must be turned so that it lies in 
the (x, z)-plane, by n'(x, y, z) = n(x,z,L — y). (Note that this only moves the vectors 
without changing their directions.) The remaining empty sub-blocks are filled with vectors 
pointing up. The above construction makes sense only if there are no discontinuities at 
the boundaries. In particular, it is necessary that also at the center of the initial un-knot 
the vector points up: n(|L, ^L, z) = (0, 0, 1). 

In Figure |3] we show the initial and final results for some linked configurations [p~2]j . On 
the left is the initial configuration (constructed in the way described) after a few iterations. 
In each case the ring on the (x, y)-plane (on the left side of the linked combination) is the 
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same, it is right-handed with Hopf charge +1 and has counterclockwise color direction 
if we look from the positive z direction. In the top two figures the part that lies in the 
(x, z)-plane also has charge +1, but the parts are linked differently, in (a) the linking is 
such that the color direction is down through the un-knot in the (x, y) plane and then the 
linking number is +2, for (b) it goes up and the linking number is —2, c.f., Figure |j. For 
the two lower figures the other un-knot is left-handed (charge —1) and when we recall the 
different definitions of direction mentioned above, the direction in (c) is from above and 
for (d) from below. 

The total charge of the linked configuration is the sum of the charges of the initial un- 
knots and ±2 for linking, and it is this total charge that is conserved under minimization. 
For (b) the sum is zero and indeed the final configuration is on its way to the vacuum 
where all vectors point up. For (c) the sum is +2 and for (d) —2 which are the charges 
of the final un-knots. It should be noted that we did not continue the minimization too 
long after the form of the final configuration became clear, because it was already known 
that deformed charge 2 un-knots return to their original shape. For (a) the sum is +4 
and the configuration seems to be stable. The total energy of this configuration is not too 
far from what is predicted by the VK bound (within 10%). A plot of energy vs. iteration 
step is given in Figure ^. 

The form of the final state depends on the relative color positions of the two un-knots. 
This is illustrated in Figure |6|. The two final configurations are different because the color 
orientations of the initial un-knots in the (x, z)-plane were different, corresponding to a 
90° rotation around the y axis [12]. The initial color difference has not changed much in 
the minimization, instead we see that the un-knots are twisted differently. 

We also studied the dynamics related to the four ways of combining un-knots of charge 
±2 [12]. The results agree with the addition rule, their energy plots have been included 
in Figure ||. Again the 2 + 2 + 2 configuration seems to be stable, when we stopped 
minimization its energy was within 7% of the VK bound. 

It should be emphasized that the minimization dynamics represents a smooth defor- 
mation of the initial configuration. One way to explain the observed deformations is to 
describe the configurations in terms of framed links and then to deform them using Kirby 
moves@, the details will be discussed elsewhere. 

Recently Battye and Sutcliffe observed that higher charge un-knots sometimes deform 
into more knotted but lower energy configurations [P!D[ . They found, e.g., that a perturbed 
charge 6 un-knot deformed into the linked combination 2 + 2 + 2 discussed above. For 
charge 4 they found a strongly deformed un-knot and it would be interesting to check 
whether its final state energy really is less than our 1 + 1 + 2 configuration. Furthermore, 
they found that the charge 7 un-knot deformed into a trefoil. In order to test whether a 
trefoil is also obtained form other initial configurations with total charge 7 we constructed 
a linked combination of charge +4 and +5 un-knots with linking number —2. The initial 
state quickly deformed into a wiggled charge 7 un-knot which then more slowly deformed 
into a trefoil. On the other hand, the combination 2 + 3 + 2 remained stable and relaxed 
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to equally low energies (within 5%). 

Our results show that the set of minimum energy configurations of Faddeev-Hopf knots 
is very rich, and the minimization can take a very complicated route. It would now be 
very interesting to characterize further the topological properties of these types of knotted 
configurations. In particular, it would be useful to know which configurations of the same 
Hopf charge can be smoothly deformed to each other, and whether in such cases there 
exists a deformation sequence where energy always decreases. In any case we know now 
that the associated knot theory must take into account the handedness and direction of 
the curve forming the knot. 
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Figure captions 



Figure 1. The final configurations for charge 1 and 2 un- knots. The displayed iso-surfaces 
correspond to the equator, i.e., have n 3 = 0. 

Figure 2. Constructing a linked combination from separate un-knots by cut and paste. 
The vector field has the vacuum direction at the boundaries marked by thick lines. 

Figure 3. Initial and final configurations for the four ways of linking a charge +1 un-knot 
with a charge ±1 un-knot. On the left-hand column we have the initial configuration, at 
the center some intermediate state, and on the right the "final" state (iso-latitude surfaces 
n 3 = —0.5). The charges and linking numbers are (a) 1 + 1 + 2, (b) 1 + 1 — 2, (c) 1 — 1 + 2, 
(d) 1 - 1 - 2. 

Figure 4. Rules for linking number from color direction. 

Figure 5. Energy as a function of the iteration step for various linked combinations. A: 
1 + 1 + 2, B: 1 + 1 - 2, C: 1 - 1 + 2, D: 1 - 1 - 2, E: 2 + 2 + 2, F: 2 + 2 - 2, G: 2 - 2 + 2, 
H: 2-2-2. 

Figure 6. Two final states for the 1 + 1 + 2 configuration. The initial states had different 
color orientation in the un-knot on the (x,z) plane (iso-latitude surfaces n 3 = —0.8). 
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Figure 1: The final configurations for charge 1 and 2 un- knots. The displayed iso-surfaces 
correspond to the equator, i.e., have n 3 = 0. 




Figure 2: Constructing a linked combination from separate un-knots by cut and paste. 
The vector field has the vacuum direction at the boundaries marked by thick lines. 
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Figure 3: Initial and final configurations for the four ways of linking a charge +1 un-knot 
with a charge ±1 un-knot. On the left-hand column we have the initial configuration, at 
the center some intermediate state, and on the right the "final" state (iso-latitude surfaces 
n 3 = —0.5). The charges and linking numbers are (a) 1 + 1 + 2, (b) 1 + 1 — 2, (c) 1 — 1 + 2, 
(d) 1 - 1 - 2. 
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Figure 4: Rules for linking number from color direction. 
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Figure 5: Energy as a function of the iteration step for various linked combinations. A: 



1 + 1 + 2, B: 1 + 1 
H: 2 - 2 - 2. 



2, C: 1 - 1 + 2, D: 1 - 1 - 2, E: 2 + 2 + 2, F: 2 + 2 - 2, G: 2 - 2 + 2, 
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Figure 6: Two final states for the 1 + 1 + 2 configuration. The initial states had different 
color orientation in the un-knot on the (x,z) plane (iso-latitude surfaces n 3 = —0.8). 
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